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THERMAL BARRIER COATING LIFE PREDICTION MODEL DEVELOPMENT* 
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The objective of this program is to establish a methodology to predict Thermal 
Barrier Coating (TBC ) life on gas turbine engine components. The approach involves 
experimental life measurement coupled with analytical modeling of relevant degrada- 
tion modes. The coating being studied is a flight qualified two layer system, de- 
signated PWA 264, consisting of a nominal ten mil layer of seven percent yttria 
partially stabilized zirconia plasma deposited over a nominal five mil layer of low 
pressure plasma deposited NiCoCrAlY . Thermal barrier coating degradation modes being 
investigated include: thermomechanical fatigue, oxidation, erosion, hot corrosion, 
and foreign object damage. 

The program is divided into two phases; the first phase, currently in progress, 
has identified cyclic mechanical damage to the ceramic and interfacial oxidation as 
the predominant degradation and failure modes (Task I) and is developing (Task II) 
and substantiating (Task III) a correlative life prediction model for these pre- 
dominant modes. Phase II will address mechanistically based modeling for all 
relevant degradation modes. 

PHASE I, TASK I - FAILURE MECHANISM DETERMINATION 

The objective of Task I was to identify predominant TBC failure modes and to 
develop a preliminary correlative life prediction model for these modes. The ap- 
proach to failure mode identification included an extensive review of experimental 
and flight service hardware, together with a laboratory test program designed to 
study the influence of driving forces such as temperature, frequency, transients, 
environment, coating thickness, etc. on degradation and failure life. 

Results of the hardware evaluation indicate the predominant failure mode to be 
thermomechanical spallation of the ceramic coating layer, resulting from the forma- 
tion of a dominant ceramic crack parallel and closely adjacent to the metal -ceramic 
interface as shown in figure 1. 

Laboratory burner rig and furnace exposure test results show cyclic coating life 
to be reduced by increased temperature, increased coating thickness, reduced cycle 
rate, pre-induced oxidation damage, and artificial sea salt injected into the burner 
flame. Analysis of these data leads to the conclusion that cyclic mechanical damage 
to the ceramic and interfacial oxidation of the underlying metal layer are the pre- 
dominant degradation modes in clean fuel cyclic thermal exposure. Exposure to salt 
contamination more severe than that normally encountered in commercial flight 
service can cause premature ceramic failure. 

Metal lographic examination of specimens removed from the burner rig test at 
various fractions of expected spallation life shows ceramic cracking as early as 20% 
of expected life. Examination of crack morphology at successively increasing life 
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fractions suggests that spallation results from progressive link-up of adjacent sub- 
critical cracks as opposed to subcritical growth of a single dominant crack as 
presented in figure 2. Direct evidence of interfacial initiated cracking is diffi- 
cult to find, indicating that the role of oxidation in damage accumulation may be 
less direct then previously thought. It is interesting to note that examples of 
scale initiated cracking are easier to find in older specimens, occurring in the 
same structure together with larger numbers of well developed longer cracks which 
appear to be isolated from the interface. This observation could suggest that the 
thicker oxide scale developed at longer exposure times can initiate cracks, but 
that this is not the "critical" damage mode in the sense that those cracks which 
propagate to failure are initiated early in life and appear to be isolated from the 
interface. An example of this is found in figure 3. 

To support the Task I preliminary modeling effort, mechanical property tests 
were conducted on bulk ceramic fabricated by plasma deposition with a structure that 
closely simulates the strain tolerant ceramic coating. The most significant result 
of these tests is the observation of highly non-linear tensile stress-strain be- 
havior at all temperatures between ambient and 2200°F. This data is summarized in 
figure 4. While there is significant variability of initial and overall stiffness, 
the basic non-linear shape of the stress-strain curve is similar at temperatures up 
to 2000°F . At 2200°F there is substantially more curvature than at lower tempera- 
tures. Both ultimate tensile strength and tensile failure strain are relatively low 
(^SKsi and « 0.3%) at all temperatures. 

Compressive stress-strain behavior, summarized in figure 5, differs signifi- 
cantly from tensile behavior; compressive strengths are much higher than tensile 
strengths, and there appears to be distinct linear and non-linear segments to the 
stress-strain curves. The 1000°F and 1600°F compressive stress-strain curves clearly 
are shaped differently than corresponding tensile stress-strain curves. At 2200°F, 
compressive deformation begins to resemble tensile deformation. Initial compressive 
stiffness appears to be essentially independent of temperature in the range studied. 

As shown in figures 6 and 7, the strain tolerant ceramic exhibits a significant 
creep response at 1800°F and 2200°F; no response was observed in a tensile creep 
test conducted at 1000°F. Based on this observation, it will be important to incor- 
porate a time dependent material response in the advanced modeling effort. 

The data plotted in figure 8 shows an apparently real fatigue response in the 
strain tolerant ceramic, but with a stress dependence substantially different from 
that observed in metals. Whereas metallic materials typically exhibit slopes ranging 
from ~ -1.5 with reversed plasticity to » -8 in the fully elastic range, the data 
in figure 8 appears to have a slope on the order of -50. Specific degradation and 
failure mechanisms responsible for this very stress sensitive fatigue behavior are 
not presently understood. 

Preliminary fracture toughness tests indicate that toughness for cracking per- 
pendicular to the splat structure is on the order of 0.5 Ksi \j in; it is expected 

that toughness for in-plane cracking, where predominant failure cracks are located 
in the coating, would be even lower. 

PHASE I, TASK I PRELIMINARY LIFE PREDICTION MODEL DEVELOPMENT 

The preliminary life prediction model focuses on the two major damage modes 
identified in the laboratory testing described above. The first of these modes in- 
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volves a mechanical driving force, resulting from cyclic strains and stresses in the 
ceramic layer caused by thermally induced and externally imposed loads. The second 
is an environmental driving force which appears, based on the experimental results, 
to be related to "oxidation damage", due to the in-service growth of a NiCoCrAlY 
oxide scale at the metal -ceramic interface. Based on the apparently "mechanical" 
mode of ceramic failure (near interfacial ceramic cracking), and on the difficulty 
in finding metal lographic evidence of a direct physical link between the growing 
oxide scale and incipient cracking in specimens exposed to a relatively small frac- 
tion of expected life, it was elected to follow the approach of Miller (ref. 1) and 
employ an existing phenomenological fatigue model (Manson - Coffin) as the basis for 
the TBC life model. In traditional form, this model relates cyclic inelastic strain 
range to a number of cycles to fatigue failure. To incorporate an environmental 
effect, the mechanical driver is analytically modified in such a way as to reduce 
the apparent fatigue strength of the ceramic layer. The use of inelastic strain 
range as a damage driver for the ceramic coating layer is considered justified in 
view of the previously mentioned nonlinearity observed in constitutive tests con- 
ducted on the strain tolerant ceramic material, including the observation of an open 
hysteresis loop in preliminary tests with reversed loading. 

The mathematical form of the model is: (Ae-j/ACf )b = Nf 


where 


Ae.j = Total cyclic inelastic strain range 

ACf = Failure strain 

Nf = Number of cycles to failure 

b = Constant 


The total cyclic inelastic strain, Aei , is the sum of the A (aAT) strain plus 
the heatup and cooldown strains, Ac c and A£^ respectively, due to the initial 
heatup and cooldown transient part of the burner rig thermal cycle: 

Ae.j = A (a AT) + Ae c + Ae H -2 (ay.s./E) 

The failure strain, Aef, is a function of the inelastic strain and is reduced 
by the strain due to the oxide thickness ratio, 6/ 6g, where 6c is the critical 
oxide thickness which will cause ceramic failure in a single thermal cycle: 

A€ f = A€ fo ( 1 - <5 / 6 c ) c + Ae i ( <5 / 6 c ) d 

The static failure strain ACf 0 , is the strain required to fail the ceramic 
in the absence of bond coat oxidation; c and d are "adjustable" constants set equal 
to unity in the preliminary analysis. Results achieved using this model to correlate 
approximately one hundred individual burner rig test results obtained on Task I are 
shown in figure 9. Achievement of a correlation coefficient of almost 0.9 is con- 
sidered quite a good fit for this preliminary attempt. 

PHASE I, TASK I - PRELIMINARY MODEL VERIFICATION TESTING 


To challenge the preliminary model, three burner rig verification tests were 
conducted using a single, internally cooled hollow specimen. This specimen permits 
exposure of the ceramic with a steady state through a thickness gradient to more 
closely simulate engine exposure of the coating, and also allows more precise in- 
strumentation and control of the thermal environment. The hollow verification test 
specimen is twice the diameter of the previously utilized specimen, and rotates 
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about its own axis to assure circumferential temperature uniformity. These substan- 
tial changes from the Task IB experimental condition assured that the preliminary 
model was effectively challenged by the verification testing. 

Comparisons between observed and predicted cyclic life, presented in figure 10 
indicate the need for significant model refinement in Task II. While the model pre- 
dicts the uncooled test result quite accurately, the two cooled test results are not 
well predicted. Refinement of the model to improve prediction capability currently 
is being performed in Task II, described below. 

PHASE I, TASK II - MAJOR MODE LIFE PREDICTION MODEL DEVELOPMENT 

The objective of this task is to improve the prediction capability of the pre- 
liminary model developed in Task I. The approach involves refinement of the analyti- 
cal model as well as the generation of additional "design data" for model calibra- 
tion, using the improved test method first employed for the Task IC verification 
tests. 

Analytical enhancements involve improved modeling of the ceramic constitutive 
and time dependent behavior, incorporation of an improved oxidation model developed 
by Dr. Robert A. Miller (which has yet to be published but which has been related 
to Pratt & Whitney through private communications with NASA), and refinement of the 
finite element calculation of temperature and stress-strain distribution. In addi- 
tion to the improved simulation of engine exposure conditions, the experimental 
program incorporates an expanded parameter envelope to cover a broader range of 
mechanical and oxidation forcing functions, as depicted in figure 11. A large por- 
tion of this test matrix has been completed and the data is being used to correlate 
a new life prediction model. 

More accurate modeling of ceramic behavior, including nonlinear stress-strain 
characteristics, assymmetric tensile and compressive response, and time dependent 
inelastic deformation, is being accomplished using a time dependent constitutive 
model developed by Walker (ref. 2). An example of the application of this model to 
the 2200°F tension and compression data is shown in figure 12. These results show 
the tensile behavior to be well modeled by the Walker approach; predictive cap- 
ability in compression is quite good up to * -1% strain; beyond this point the pre- 
diction departs significantly from observed behavior. An example of a hysteresis 
loop calculated for a typical strain emphasis burner rig cycle is shown schemati- 
cally in figure 12- There are two significant observations concerning the calculated 
behavior. First, compressive strain remains below the -1% limit, beyond which the 
prediction capability of the model breaks down. Second, the model predicts a very 
open loop with a quite large reversed plastic strain range to drive ceramic fatigue 
damage. Prediction for all of the experimental burner rig cycles currently is being 
made with this model for incorporation into the improved life prediction system. 

Oxidation experiments were conducted using the program's substrate - TBC system, at 
the NASA Lewis Research Center. Dr. Robert Miller used this data to develop a new 
more accurate oxidation model for the improved life prediction system. The recom- 
mended weight gain expression for the data showed that the oxide growth rate is not 
parabolic, but approximates t 0 - 29 . This new oxide growth rate expression in 
terms of oxide thickness is being incorporated into the advanced life prediction 
model. 
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PHASE I - TASK III - MODEL VERIFICATION 


The objective of this task is to validate the model developed in Task II through 
a series of approved benchmark engine mission simulation tests. Experimental ap- 
proaches to further improve the simulation capability of the test method (i.e., 
generate higher heat fluxes and larger through thickness gradients), and thus more 
effectively challenge the model, currently are being explored for incorporation in 
the Task III verification testing. Based on results of these tests, recommendations 
for further research or refinement required to arrive at a fully satisfactory engine 
life prediction methodology shall be made, if necessary. 
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Figure 1 Typical Thermal Barrier Coating Engine Failure Mode 
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Figure 2 Thermal Barrier Coating Damage Progression 


ORIGINAL PAGE IS 

OF POOR QUALITY 




476 



Figure 2 (Continued) 


477 



1 00% (180 hrs) 


Figure 2 (Continued) 
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Figure 3 BSI of TBC After 105 Hours of 
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Rig Test Time at 2100°F/Short Cycle/Fast Heatup 





ENGINEERING TENSILE STRAIN, PERCENT 


Figure 4 Representative Strain Tolerant Ceramic Tensile Stress Strain Curves at 
Various Temperatures. Room temperature strain data measured by strain 
gage; temperature curves obtained from corrected cross head displacement. 



ENGINEERING COMPRESSION STRAIN, PERCENT 


Figure 5 Representative Strain Tolerant Ceramic Compressive Stress-Strain Curves 
at Various Temperatures. Compressive strains calculated from corrected 
crosshead displacement. 
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Figure 8 S-N Curve for 7YsZ (1000°F and 1600°F data plotted together) 
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Figure 9 Preliminary Life Model Correlation 
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Figure 12 Compression and Tension Data at 2200°F Kevin Walker Model Predictions 
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Figure 13 


Schematic of Thermally Driven Ceramic Stress 
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